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1. Introduction 

We have developed a technique for constructing addition theorems for the azimuthal 
Fourier coefficients of fundamental solutions for linear homogenous partial differential 
equations on d- dimensional Riemannian manifolds. In our technique, one must 
construct azimuthal Fourier expansions for this fundamental solution and compare it 
with eigenfunction expansions for that fundamental solution in particular separable 
rotationally-invariant coordinate systems. The construction of eigenfunction expansions 
for fundamental solutions in separable coordinate systems is in general non-trivial. 

In this paper, we apply this technique to generate addition theorems from 
eigenfunction expansions for a fundamental solution of the polyharmonic operator on 
(i-dimensional Euclidean space R d in Vilenkin-Kuznetsov-Smorodinskh polyspherical 
coordinate systems [251 126] (hereafter Vilenkin) . These coordinate systems are subgroup 
type coordinate systems, as they can be descibed using a maximal subgroup chain, and 
are perhaps the most studied rotationally-invariant coordinate systems in existence for 
this operator. Here we have computed azimuthal Fourier expansions for a fundamental 
solution of this operator, as well as the corresponding eigenfunction expansions in 
polyspherical coordinate systems. In each case, the comparison of these two expansions 
yields new addition theorems for the azimuthal Fourier coefficients. 

A fundamental solution of the polyharmonic equation on R d is given either in 
terms of solutions which have power-law or logarithmic behavior as a function of the 
distance between two points on this space (see for instance Schwartz (1950) ([221 P- 45]); 
Gel'fand & Shilov (1964) ([TJ], p. 202)). In a recent paper (Cohl & Dominici (2010) 
[5]), we derived a complex identity which determined the Fourier coefficients of a 
power-law fundamental solution of the polyharmonic equation. The Fourier coefficients 
were seen to be given in terms of associated Legendre functions. The eigenfunction 
expansions for a power-law fundamental solution of the polyharmonic equation in 
a particular polyspherical coordinate system can be derived using a Gegenbauer 
polynomial expansion for the relevant kernel (see Corollary 15.41 below) in conjunction 
with the addition theorem for hyperspherical harmonics. This addition theorem is given 
by 

C d J 2 -\co^) = 2(d ~ 2 \ 7r ^ 2 |/ y Y«(5i)Y?W), (1) 
v (2n + d - 2)r(d/2) ^ v ' v ; v ; 

(for a proof see Wen & Avery (1985) [27]; see also section 10.2.1 in Fano & Rau (1996) 
[11] ) where K stands for a set of (d — 2)-quantum numbers identifying degenerate 
harmonics for a given value of n G N , and cos 7 is the cosine of the separation angle 
between two arbitrary vectors x, x' G R d (see ( I12j) below). The functions : S^" 1 — > C 
are the normalized hyperspherical harmonics, with x, x' G S d_1 being unit vectors in 
the direction of x, x' G R d respectively. The Riemannian manifold is defined as the 
set of all points in R d+1 such that Xq + . . . + x\ = R 2 {R > 0), with the metric induced 
from that of the ambient Euclidean space. We denote the rf-dimensional hypersphere of 
unit radius by S d := Sf. 
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There will be a separate addition theorem corresponding to each rotationally- 
invariant coordinate system. It should be mentioned that the large number and 
complexity of coordinate systems increases multi-fold when one increases dimension. 
One could also allow for the possibility of non-subgroup type rotationally-invariant 
coordinate systems such as parabolic and spheroidal coodinates but these will not 
be addressed in the present paper. See for instance Kalnins (1986) [T7] and Cohl et 
al. (2000) [8] for some examples of this type of coordinate systems connected with the 
Laplace operator on R 3 . In this paper we give examples of addition theorems from a 
power-law fundamental solution of the polyharmonic operator in Euclidean space using 
Vilenkin's polyspherical coordinate systems in low dimensions (i.e., d G {3,4}) and also 
in higher dimensions. 

This paper is organized as follows. In section [2] we describe the kernels associated 
with a fundamental solution of the polyharmonic equation on Euclidean space R d . 
In section [3] we describe Vilenkin's polyspherical coordinate systems. In section 
H] we describe the normalized hyperspherical harmonics associated with Vilenkin's 
polyspherical coordinates. In section [5] we prove several theorems associated with Jacobi, 
Gegenbauer, and Chebshev polynomial expansions of the kernels associated with power- 
law fundamental solutions of the polyharmonic equation on Euclidean space. In section 
[6] we derive and discuss addition theorems in Vilenkin's polyspherical coordinates for 
the Fourier coefficients of a fundamental solution for the polyharmonic equation on R d . 

Throughout this paper we rely on the following definitions. Let a 1 ,a 2 ,a 3 , . . . G C, 
with C being the set of complex numbers. If i,j G Z and j < i, then YH > n =% a n = 
and Yln=i a n — 1- The set of natural numbers is given by N := {1,2,3, . . .}, the set 
N := {0, 1, 2, . . .} = N U {0}, and Z := {0, ±1, ±2, . . .}. The sets Q and R represents 
the rational and real numbers respectively. 

2. Fundamental solution of the polyharmonic equation and the 
non-logarithmic Fourier series 

If $ : R d — y R satisfies the polyharmonic equation given by 



where x G R d , A : C p (R d ) — > C p 2 (R d ) for p > 2, is the Laplacian operator defined by 



where we take p to be an integrable function so that a solution to fl3]) exists. A 




(2) 
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where 5 is the Dirac delta function (generalized function) and x' G R d . The Euclidean 
inner product (-, •) : R d x R d — > R defined by 

(5) 



R d -> [0, oo), on the finite-dimensional 



(x, x ) :— + . . . + x d x dl 

induces a norm (the Euclidean norm) || ■ 
vector space R d , given by ||x|| := y/ (x, x). 

In the rest of this paper, we will use the gamma function T : C \ — N — > C, which 
is a natural generalization of the factorial function (see for instance Chapter 5 in Olver 
et al. (2010) [5T] ). A fundamental solution of the polyharmonic equation is given by the 
following theorem. 

Theorem 2.1. Letd,keN. Define Gt:(R d x R rf ) \ { (x, x) : x G R d } — > R such that 



£fc(x,x) := < 



>fc+d/2+l 



/||2fc-d 



(jfe- 1)! (A;-d/2)!2 2fc " 1 7r cZ/2 



T(d/2- fc)||x-x 



l\\2k-d 



(k - 1)! 2 2k TT d/2 
where (3 p ^ d £ Q is defined as (3 p ^ d := | [i? p 



(log||x-x'|| - P k - d /2,d) 

if <i even, /c > c//2, 
otherwise, 

H d /2+ P -i — H d / 2 -i] , Hj being the jth 



harmonic number 



j 1 



i=l 



then Q d is a normalized fundamental solution for (— A) k on R d . 

Proof. See Cohl (2010) [4] and Boyling (1996) [3]. 

We consider parametrizations of points on R d which are given in terms of coordinate 
systems whose coordinates are curvilinear, i.e., based on some transformation which 
converts the Cartesian coordinates to a coordinate system with the same number of 
coordinates in which the coordinate lines are curved. We consider solutions of the 
polyharmonic equation ([2]), in a curvilinear coordinate system, which arises through the 
theory of separation of variables. We refer to coordinate systems which yield solutions 
through the separation of variables method as separable. In this paper we restrict our 
attention to separable rotationally-invariant coordinate systems for the polyharmonic 
equation on R d which are given by 

Xi = ^i,---,^-i)cos0 

%2 = #(£i,---,£d-i)sin0 

^3 = X3(€l,'--,€dr-l) >• (6) 



x d — x d\£,li ■ ■ ■ > Cd-l) 

These coordinate systems are described by d coordinates: an angle G [0, 2ir) plus 
(d — 1) -curvilinear coordinates (£ 1; . . . , Rotationally-invariant coordinate systems 
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parametrize points on the (d — l)-dimensional half-hyperplane given by = const 
and R > using the curvilinear coordinates • • • , £<z_i)- A separable rotationally- 
invariant coordinate system transforms the polyharmonic equation into a set of d- 
uncoupled ordinary differential equations with separation constants m G Z and kj G R 
for 1 < j < d — 2. For a separable rotationally-invariant coordinate system, this 
uncoupling is accomplished, in general, by assuming a solution to (j2J) of the form 

d-l 

= e w ^(Ci,...,^_i) Aifam, k t , . . . , k d ^ 

=i 



where the properties of the functions 1Z and A t , for 1 < i < d — 1, and the constants 
fcj for 1 < j < d — 2, depend on the specific separable rotationally-invariant coordinate 
system in question. Separable coordinate systems are divided into two separate classes, 
those which are simply separable (1Z = const), and those which are ^-separable. For a 
general description of the theory of separation of variables see Miller (1977) [20] . 

The Euclidean distance between two points x, x' G R d , expressed in a rotationally- 
invariant coordinate system, is given by 

^[x-cos^-^)] 172 , (7) 



x — X 



where the toroidal parameter \ > 1, is given by 



d 

)2 , r>/2 , V~Y _ '\2 



R 2 + R' 2 + J2( 



X '' 2RR' ' (8) 

where R,R' G (0, oo) are defined in (jSJ) for x, x' G R d . The hypersurfaces given by 
X > 1 equals constant are independent of coordinate system and represent hypertori of 
revolution. 

From Theorem 12.11 we see that, apart from multiplicative constants, the algebraic 
expression l k : (R d x R d ) \ {(x, x) : x G R d } — > R of an unnormalized fundamental 
solution for the polyharmonic equation on R d for d even, k > d/2, is given by 

lf(x, x') := ||x - x'|| 2fc " d (log ||x - x'|| - k _ d/2td ) . (9) 

By expressing if in a rotationally-invariant coordinate system OH]) we obtain 



^log (2RR')-/3 Pil 



[x ~ cos(0 - < 

(2RR') P [ X - cos(0 - 0')] P log [X ~ cos(0 - 0')] , (10) 



^(x,x')= (2RR 1 ) 1 
1 

+ 2 

where p = k — d/2 G N . For the polyharmonic equation on even- dimensional Euclidean 
space R d with 1 < k < d/2 — 1, apart from multiplicative constants, the algebraic 
expression for an unnormalized fundamental solution of the polyharmonic equation 
hf : (R d x ~R d ) \ {(x, x) : x G R d } -> R is given by 

r di i\ ii i\\2k—d 

f)fe(x,x) := x-x 
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By expressing fyf in a rotationally-invariant coordinate system we obtain 

f^(x,x') = (2RR!y q [ X - C os{<P-<p')Y\ (11) 

where q = 2k — d G N. 

By examining ([TO]) . ffTTj) . we see that for computation of Fourier expansions about 
the azimuthal separation angle (0 — <p') of if and t)f, all that is required is to compute 
the Fourier cosine series for the following three functions f x , h x : R — > (0, oo) and 
g x : R — > R defined as 

/ X W : = (X - cos^f, 

9 X ($) '■= (X~ cos ^) p log (x - cos ^) , 

fyeWO := (* ^ cosV')" 9 , 
where p G N , q G N and x > 1 is a fixed parameter. In a polyspherical coordinate 
system, the separation angle 7 G [0,7r], is the smallest angle measured between two 
arbitrary non-zero length position vectors x, x' G R fi . It is defined through the relation 

ll x llll x II 

using the Euclidean inner product and norm (cf. (JSJ)). 

The Fourier series of f x is given in Cohl & Dominici (2010) [5] (cf. (4.4) therein), 
namely for p G N , 

(, - x y = ( 2 2 - if' 2 £ £ ° ( "^)7"" p p (^f=y) T »<-). < 13 > 

where e n G {1,2} is the Neumann factor defined by e n := 2 — <5 n , 5 n G {0, 1} is the 
Kronecker delta, and the Pochhammer symbol (rising factorial) (•)„ : C — > C for n G N 
is defined by 



[ Z )n 

i=l 

We have used Whipple's formula in f fT3"j) (see for instance Olver et al. (2010) J2TJ 
(14.9.16), (14.9.17)]) to convert the associated Legendre function of the second kind Q„ 
appearing in [5] to the associated Legendre function of the first kind P£. The associated 
Legendre function of the second kind Q% : C\ (— 00, 1] — > C, u + fi — N, can be defined 
in terms of the Gauss hypergeometric function as follows (Olver et al. (2010) J2TJ, (14.3.7) 
and section 14.21]) 

for |z| > 1 and by analytic continuation of the Gauss hypergeometric function elsewhere. 
Similarly, the associated Legendre function of the first kind can be defined using the 
Gauss hypergeometric function ((14.3.6) and §14.21 (i) in Olver et al. (2010) [21] ) 

* W:= * (i±lY'\ Fl (-^^^), (14) 
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where |1 — z\ < 2, and elsewhere by analytic continuation. The Gauss hypergeometric 
function 2 i*i : C x C x (C \ — N ) x {z G C : \z\ < 1} — > C can be defined in terms of 
the following infinite series 

( a,b \ ^ (a) n (6) n z n 

(see for instance Chapter 15 in Olver et al. (2010) [21]), and elsewhere by analytic 
continuation. 

The Fourier series of h x is given in Cohl & Dominici (2010) [5] (4.5)] (we have used 
the Whipple formula), namely for p G N, 

i ( 7 2 _ ir p/2 00 ( 7 \ 

i^W =W^'" +p - 1)!P - (t^i) W (16) 

In order to compute Fourier expansion of if ffTUj) in separable rotationally-invariant 
coordinate systems, all that remains is to determine the Fourier series of g x (see Cohl 
(2012) [3]). We will leave the discussion of Gauss hypergeometric orthogonal polynomial 
expansions for logarithmic fundamental solutions of the polyharmonic equation on R d 
(from g x ) to a later paper and from here on out we will only discuss expansions of 
power-law fundamental solutions (from f x , hJ) of the polyharmonic equation. 

3. Vilenkin's polyspherical coordinate systems on H d 

In this section we review results in the literature concerning Vilenkin's polyspherical 
coordinates, which are subgroup type coordinates. Subgroup type coordinate systems 
on H d are coordinate systems which can be described by a maximal chain of subgroups 
of the Euclidean group E(d). 

Polyspherical coordinates are specific subgroup type coordinate systems which 
are described by a radial coordinate r G [0, oo) plus (d — l)-angles which together 
parametrize points on a (d — 1) -dimensional hypersphere. First we will discuss a 
general procedure for constructing polyspherical coordinate systems called the "method 
of trees," and then give some examples of polyspherical coordinates which will be used 
in the rest of the paper. 



3.1. Method of trees 

Subgroup type coordinates are coordinates in which the very useful method of trees 
(a graphical method for generating coordinate transformations for subgroup type 
coordinates) can be adopted. Corresponding to each of these subgroup chains are 
a set of tree diagrams. Each of these trees corresponds to a particular separable 
coordinate system in <i-dimensions. See Izmest'ev et al. (1999) [15] (and below) for 
explicit parametrizations for the corresponding coordinate systems and trees. Describing 
polyspherical coordinate systems in terms of rooted trees was originally developed by 
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Vilenkin, Kuznetsov & Smorodinskii (1965) (see also Vilenkin (1968) (25]) and has 
since been used extensively by others in a variety of contexts (see for instance Izmest'ev 
et at (1999) p], Izmest'ev et al (2001) p], Kil'dyushov (1972) [18]). 

In these rooted trees, there are two types of nodes, leaf nodes and branching nodes. 
For a coordinate system on R d , there are d leaf nodes, each corresponding to the 
particular Cartesian component of an arbitrary position vector x £ R d . The branching 
nodes split into two separate branches, one up to the left and one up to the right. Each 
branch emanating from a branching node will end on either a leaf node or on another 
branching node. There are four possibilities for branching nodes (see figured]). 




Figure 1. This figure shows the possibilities from left to right for branching nodes 
of type a, b, b' and c. For type a, both branches end on a leaf node. The angle 
corresponding to this type of branching node is <fi a £ [0, 27r). For type b, the left 
branch ends on a leaf node and the right branch ends on a branching node. The 
angle corresponding to this type of branching node is 9^ £ [0, n]. For type b', the left- 
branch ends on a branching node and the right branch ends on a leaf node. The angle 
corresponding to this type of branching node is 6 b > £ [—w/2, n/2]. For type c, both the 
left and right branches ends on branching nodes (cells of type c are only possible for 
d > 4). The angle corresponding to this type of branching node is d c <G [0, ir/2]. 

A linear partial differential equation on R d may admit solutions via the separation 
of variables method. If the method works, then it produces d-ordinary differential 
equations with (d — l)-separation constants which are determined by the conditions 
imposed on the problem being solved. The sets on which the separation constants are 
contained depend on the geometric structure of the problem being solved. For instance, 
separation of variables in polyspherical coordinates with (d — l)-angles, for Laplace's 
equation on R d produces (d — l)-separation constants, each of which are called quantum 
numbers. The quantum numbers corresponding to these angles are all integers. This is 
due to the fact that a hypersphere is a compact manifold. If the separation constants 
are integers, then they are referred to as quantum numbers. Quantum numbers for 
a particular tree label the basis of separable solutions for Laplace's equation in that 
particular coordinate system. With each branching node of the tree, we associate a 
quantum number as well as an angle. 

The quantum number corresponding to a (27r)-periodic (azimuthal) angle is called 
an azimuthal quantum number. Each azimuthal angle corresponds to a branching node 
of type a, and an azimuthal quantum number m £ Z. A natural consequence of the 
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method of trees is that there must exist at least one azimuthal angle for each tree, and 
therefore also for each polyspherical coordinate system. 

Branching nodes of type b, b' and c, as described above, are associated with 
angles, which in turn are associated with quantum numbers which we refer to as angular 
momentum quantum numbers I G N . The naming of these quantum numbers is due 
to the intimate connection between the properties of the rotation group SO(d) and 
the quantum/classical description of angular momentum in physics. See for instance 
Chapter 10 in Fano & Rau (1996) [II]. 

There is always at least one branching node, the root, and all branching nodes 
correspond to a particular angle and quantum number. Let us associate each branching 
node with an angle and its corresponding quantum number. Parametrization of 
polyspherical coordinate systems are accomplished as follows. Starting at the root, 
traverse the tree upward until you reach the leaf node corresponding to x { . The 
parametrization for is given by the hyperspherical radius r multiplied by cosine or sine 
of each angle encountered as you traverse the tree upward until you reach the leaf node 
corresponding to Xj. If you branch upwards to the left or upwards to the right at each 
branching node, multiply by the cosine or sine of the corresponding angle respectively. 
This procedure produces the appropriate transformation from polyspherical coordinates 
to Cartesian coordinates. 

There are large numbers of equivalent trees and an even larger number of possible 
trees, each with their own specific polyspherical coordinate system. The enumeration 
of these trees are characterized as follows. For d, b d G N, let b d be the total number of 
trees (polyspherical coordinate systems) on R d . Then b 1 = 1 is the number of possible 
1-branch trees. A 1-branch tree corresponding to a polyspherical coordinate system 
does not exist in isolation. However, the concept of a 1-branch tree is useful. For 
instance, one can generate a new tree from a pre-existing tree by adding a single branch 
to either the left or to the right of a pre-existing tree and in order to perform the correct 
enumeration, one must incorporate 1-branch trees. The following recurrence relation 
gives the total number of possible polyspherical coordinate systems (trees) for arbitrary 
dimension 



Using the recurrence relation given by ( fl7l) . the first few elements of the sequence are 
given by 

(b d : d G {2, . . . , 13}) = (1, 2, 5, 14, 42, 132, 429, 1430, 4862, 16796, 58786, 208012). 

These numbers, or the way to count these types of trees, are given in terms of the 
Catalan numbers C n (see for instance Sloane Integer Sequence A000108 [23] or p. 200 
in Stanley (1999) [21]), i.e., 



d-l 




(17) 



i=l 
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where 



n 



nl 



kj k\(n-k)\ 

is the binomial coefficient for k,n e N with < k < n. 

If a d e N is the total number of equivalence classes for equivalent trees for each 
dimension d, then the following recurrence relation gives the number of equivalence 
classes for arbitrary d 



Yd/2\ 

E 

i=l 



a i a d- 



d/2-1 

aiCLd - 



+ ^d/2 {ad/2 + 1) 



if d odd 



if d even. 



Using the recurrence relation given by 
given by 



181) . the first few elements of the sequence are 



(a d : d e {2, . . . , 13}) = (1, 1, 2, 3, 6, 11, 23, 46, 98, 207, 451, 983). 

These numbers, or the way to count these types of trees, are given in terms of the 
Wedderburn-Etherington numbers (see for instance Sloane Integer Sequence A001190 
[23J). The equivalence classes are established by a left-right symmetry in the topology 
of the trees. The types of polyspherical coordinates that one can construct, for a given 
dimension, are described in Izmest'ev et al. (1999, 2001) [T5l IT6] . 

We use a left-to-right recursive naming language for our trees based on a depth-first 
search (see pp. 540-549 in Cormen et al. (2001) [9]). This naming language is given by 
listing the types of branching nodes available in a particular tree. 



3.2. Vilenkin's polyspherical coordinate systems 

A particular polyspherical coordinate system partitions R d into a family of concentric 
(d — 1) -dimensional hyperspheres, each with a radius r G (0, oo), on which all 
polyspherical coordinate systems for Sf" 1 may be used. One must also consider the 
limiting case for r = to fill out all of R d . The hyperspherical radius is computed using 
Cartesian coordinates through 




i=i 



A polyspherical coordinate system yields solutions to Laplace's equation on R through 
the separation of variables method. 

The method of trees can be used to construct the cosine of the separation angle 
ffT2"]) in a direct manner. The cosine of the separation angle will be given by the sum 
of rf-terms, each corresponding a leaf node of the tree. There is a unique path starting 
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from the root to each leaf node. The cosine of the separation angle can be constructed 
using the following formula 

d Ni 

cos7 = J^n^-WOAiWi)' 

i=l 3=1 

where N { is the number of branching nodes encountered from the root to the leaf 
node, A±j : R — > [—1,1] is either the trigonometric cosine or sine function depending 
respectively on whether the left branch or right branch is chosen, and i/)^ G R is the 
angle corresponding to the jth branching node for each ith leaf node. The formula for 
the cosine of the separation angle is unique for each tree. 

3.3. Low dimensional examples of Vilenkin's poly spherical coordinate systems 

The simplest example of a polyspherical coordinate system on R d occurs for d = 2 (polar 
coordinates) where there is one branching node (the root) and two leaf nodes (see figures 
[2k andEK). The left-branch ends on the leaf node corresponding to x 1 and the right 
branch ends on the leaf node corresponding to x 2 . Therefore for d = 2 we have the 
following coordinate system whose transformation formulae to Cartesian coordinates is 
given by 



x 1 = r cos 
x 2 = r sin </> 




where r G [0, oo) and <fi G [0,27r). This coordinate system is the only member of its 
equivalence class (cf. ( ITS]) ). Using our naming procedure based on a depth-first search, 
we refer to this tree as type a. These coordinates are adapted to the canonical subgroup 
chain E(2) D 0(2). The cosine of the separation angle is given by 

cos 7 = cos(0 — 0'), (20) 

and corresponding to the angle <fi is the quantum number m G Z (see figured). 



. 0C - jC OC 0C 0C jC oc 

1 2 1 2 3 1 2 : 




(a) (b) (c) 

Figure 2. Tree diagrams for two and three dimensional polyspherical coordinate 
systems of type: (a) a, (b) ba, (c) b'a. 

In d = 3 there are 2 possible topological trees, each corresponding to one of 
two different trees. Both coordinate systems are adapted to the canonical subgroup 
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chain E(3) D 0(3) D 0(2). The first tree (see figure [2)3 ) corresponds to the following 
coordinate system whose transformation formulae to Cartesian coordinates are given by 

X\ — r cos 9 1 

x 2 = rsin#cos0 > , (21) 
x 3 = rsin#sin0 J 

where 9 G [0,7r], e [0,27r). We call this polyspherical coordinate system standard 
spherical coordinates. Using our naming procedure based on a depth-first search, we 
refer to this tree as type ba. The cosine of the separation angle is given by 

cos 7 = cos 9 cos 9' + sin 9 sin 9 cos (0 — 0'), (22) 

and corresponding to the angles 9 G [0, ix] and G [0, 2n) are the quantum numbers 
/ G N and m G Z respectively 

The second tree (see figure [2b) corresponds to the following coordinate system whose 
transformation formulae to Cartesian coordinates are given by 

Xi — r cos 9 cos "| 

x 2 = r cos 9 sin > , (23) 
x 3 = r sin J 

where 9 G [— 7r/2, 7r/2], G [0, 27r). This tree is of type b'a. The cosine of the separation 
angle is given by 

cos 7 = sin 9 sin 9' + cos 9 cos 9' cos (0 — 0'), (24) 

and corresponding to the angles 9 and are the quantum numbers / G N and m G Z 
respectively. The coordinate systems of type ba and b a belong to the same equivalence 
class. 

In d — 4 there are five possible topological trees for polyspherical coordinates, each 
corresponding to a different tree. The first four coordinate systems are adapted to the 
canonical subgroup chain E(A) D 0(4) D 0(3) D 0(2). The first tree (see figures [3^ 
and H|) corresponds to the following coordinate system, whose transformation formulae 
to Cartesian coordinates are given by 

x 1 = rcos9i 

x 2 = r sin 9 1 cos 9 2 

x 3 = r sin 9 1 sin 9 2 cos 

x A = r sin 9 1 sin 9 2 sin 

where 9i,9 2 G [0,7r]. Using our naming procedure based on a depth-first search, this 
tree is of type b 2 a. The cosine of the separation angle is given by 

cos 7 = cos 9 1 cos 9[ + sin 9 1 sin 9[ (cos 9 2 cos 9' 2 + sin 9 2 sin # 2 cos(0 — 0')) . (26) 



(25) 
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] [ 





(a) 



(b) 



(c) 



OC-^ OC^ OC/^ 



0C-\ OC^ oc^ 





(d) 



(e) 



Figure 3. Tree diagrams for four dimensional polyspherical coordinate systems of 
type: (a) b a, (b) bb'a, (c) b'ba, (d) b' a, (e) ca . 



The second tree (see figure [3b) corresponds to the following coordinate system 
transformation formulae to Cartesian coordinates are given by 

Xi — r cos 9 l 

x 2 = r sin 9 X cos 9 2 cos 

x 3 = r sin 9 X cos 9 2 sin </> 

x 4 = r sin 9\ sin 9 2 

where 9 1 G [0,7r] and 9 2 G [— 7r/2, 7r/2]. This tree is of type bb'a. The cosine 
separation angle is given by 

cos 7 = cos 9 1 cos 9'i + sin 9 X sin 9[ (sin #2 sin #2 + cos 9 2 cos # 2 cos(0 — (/>')) ■ 

The third tree (see figure [3b) corresponds to the following coordinate system, 
transformation formulae to Cartesian coordinates are given by 

x 1 = r cos 9 1 cos 9 2 

x 2 = r cos 9 1 sin 9 2 cos 

x 3 = r cos # x sin # 2 sin 

x 4 = r sin 

where 9 1 G [— 7r/2,7r/2] and # 2 £ [0, 7r]. This tree is of type b ba. The cosine 
separation angle is given by 

cos 7 = sin 9 1 sin 9[ + cos 6^ cos 9[ (cos 6* 2 cos #2 + si n #2 sm $2 cos(0 — 4 > ')) ■ 



whose 



(27) 

of the 

(28) 
whose 



(29) 

of the 
(30) 
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(31) 



The fourth tree (see figure 01) corresponds to the following coordinate system, whose 
transformation formulae to Cartesian coordinates are given by 

X\ — r cos 9 1 cos 2 cos 

x 2 = r cos 9 1 cos 9 2 sin 

x 3 = r cos 9 1 sin 9 2 

x A = r sin 9 2 

where 1; 9 2 G [— tt/2, tt/2]. This tree is of type b /2 a. The cosine of the separation angle 
is given by 

cos 7 = sin 9 1 sin 9[ + cos 9 1 cos 9[ (sin 9 2 sin 9' 2 + cos 2 cos $2 cos(0 — 0')) . (32) 
The fifth tree (see figures [3b and \5jp) corresponds to the following coordinate system 
(Hopf coordinates) whose transformation formulae to Cartesian coordinates are given 
by 



•''1 

X 1 



r cos i? cos 
r cos •§ sin X 
r sin d cos 2 
r sin •& sin 2 



(33) 



where i9 G [0, 7r/2] and 4>i,4> 2 G [0, 27r). This tree is of type ca 2 . This coordinate system 
is adapted to the canonical subgroup chain £"(4) D 0(4) D 0(2) ® 0(2). The cosine of 
the separation angle is given by 

cos $ cos cos^x — 0i ) + sin ^ sin $ cos(0 2 — 02)- (34) 



cos 7 



5.^. High dimensional examples of Vilenkin's poly spherical coordinate systems 

There are many choices for polyspherical coordinates on R d , suitably defined for any 
number of dimensions d > 2. Note that these also give polyspherical coordinate systems 
on S d_1 when r = 1. The simplest example of a polyspherical coordinate system 
on R d , valid for high dimensions, is what we will refer to as standard polyspherical 
coordinates (see figure H]). Standard polyspherical coordinates are a generalization 
of (12T1 type ba) and (1251 type b 2 a). This coordinate system (and all those which 
belong to its same equivalence class) is adapted to the canonical subgroup chain 
E(d) D 0(d) D 0(d — 1) D ■ ■ • D 0(2). The transformation formulae to Cartesian 
coordinates for standard polyspherical coordinates are given by 
x 1 — r cos 9 1 
x 2 = r sin 9 1 cos 9 2 

r sin 9 1 sin 9 2 cos 9 3 



•>'3 



Xd-2 = r sin 0x ■ • • sin d _ 3 cos 6» d _ 2 

= r sin X • ■ • sin d _ 3 sin rf _ 2 cos 1 
= r sin X • • • sin 9 d _ 3 sin 9 d _ 2 sin c 



(35) 
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where 9 { G [0, vr] for i G {1, . . . , d — 2} and G [0, 27r). Using our naming procedure 
based on a depth-first search, this tree is of type b d ~ 2 a. These are the standard 
polyspherical coordinates that are usually adopted in mult i- variable calculus and in 
physical applications (see for instance Fano & Rau (1996) [H]). In these coordinates 
the cosine of the separation angle is given by 

d-2 i-1 d-2 

cos 7 = cos #jCos 9'i Y\ sin fysin 9j + cos(</> — <//) Y\ srn #jSin Q\. 

i=l j=l i=l 




Figure 4. Tree diagram for cZ-dimcnsional standard polyspherical coordinates of type 
b _2 a (standard polyspherical coordinates). 

Another example of a polyspherical coordinate system which is valid for large 
dimensions is what we will refer to as generalized Hopf coordinates (see figure |5]). These 
coordinates, valid on R 2 for q > 1, generalize two-dimensional polar coordinates ( 1T91 
type a) and four-dimensional Hopf coordinates ( 133| type ca 2 ). These coordinates are 
unique in that they correspond to the only trees which contain only themselves in 
their equivalence class (see ( |T8l) ). These coordinate systems have separated harmonic 
eigenf unctions which are given in terms of complex exponentials, and for q > 2, non- 
symmetric Jacobi polynomials (see figure [5]). These coordinates are adapted to the 
canonical subgroup chain 

E(d) D 0(d) d 0(d - 1) ® 0(d - 1) 



D 0(d - 2) ® 0(d - 2) (2) 0(d - 2) (8) 0(d - 2) D • • • D 0(2) ® • • • <g> 0(2). 
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x. 



X. 



Xy y y y y y y 
^ "3 "4 "g ~^-g "-c 




(a) 



(b) 



(c) 



Jy y y y y y y y y y y y y y y 

^ i/V^ § **^5 "'^Q *^^9 xA'-pi q »/V^ -j 2 lA^pj 3 •i'V^ ^ t/V^ ^ g 

l>1 y >. ^2 , v <l>3 , v 4>4 „ K 4>5 ^ „ <f»£y ^4>7 




(d) 

Figure 5. This figure is a tree diagram for polysphcrical generalized Hopf coordinates 
of type V 2 i on R with q = 1,2,3, 4 for (a), (b), (c), and (d) respectively. The first 
(2 9_1 — l)-branching nodes are of type c which correspond to the angles $j € [0, f ] 
and quantum numbers 1+ € N . The following (2 9 )-branching nodes are of type a 
which correspond to the angles ^ £ [0, 2ir) and quantum numbers m, £ Z. These 
coordinates correspond to Transformation Q36p . 



This coordinate system is suitably defined for dimensions d 
transformation formulae to Cartesian coordinates are given by 



%2 



X^q — I 



X„q-1 



X i 



COS r &„q-2 COS 0! 



cos $,-,9-2 sin 0! 



r COS COS $ 2 cos ^4 cos $8 
r COS COS $ 2 cos ^4 COS ^8 

r cos di sin $ 2 sin $ 5 sin $ n 
r cos i?! sin d 2 sin i? 5 sin $ n 
r sin $ x cos t9 3 cos t? 6 cos i? 12 • • • cos ^ 3 . 2 i-s cos 
r sin $ x cos t? 3 cos t? 6 cos $ 12 • • • cos $ 3 . 2 q~3 sin 



■ sin ^3. 2 9-3_ 1 cos (f> 2 q-2 



• sin-i? 



3-2 



q-3 



-l sm r 2 



9-2 



9-2 



+1 



2^-1 



r sin ^ sin $ 3 sin $ 7 sin $ 17 ■ • • sin $ 9 9-i , cos </> 9 <j-i 



2 9 for g > 1. The 



(36) 
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where ^ 6 [0, §] for % = 1, . . . , 2 q ~ 1 - 1 and & G [0, 2tt) for i = 1, . . . , 2 9_1 . Using our 
naming procedure based on a depth-first search (see section [3TTT) . these coordinates are 
of type 

V 2 9 = cV 2 9-iV 2 «-i, 

where q G N, with type V 2 = a (polar) coordinates on R 2 (see figures [2k. and[5h). The 
first non-trivial example of generalized Hopf coordinates is V 4 = cV 2 V 2 = ca 2 (Hopf 
coordinates, see figures [3b and [5b). Other examples of generalized Hopf coordinates 
are of type V 8 = cV 4 V 4 = cca 2 ca 2 and of type V 16 = cV 8 V 8 = ccca 2 ca 2 cca 2 ca 2 
coordinates (for corresponding tree diagrams, see figures [5b and EH). 

The cosine of the separation angle in this coordinate system may be given as follows. 
Consider the symbol q G r s G [—1,1] with < s < q and 1 < r < 2 q — 1. If this symbol is 
given recursively by 

q G r s = cos$ r _ 1+29 - s cos 0(._ 1+2 «-* gG 2 ^ 1 + sintf._ 1+29 - s sin#_ 1+29 _ s q Q 2 JLi, 
with q G l = 1, then the cosine of the separation angle is given by 
cos 7 = q G l q 

(cf. ( l20l . ( l34l) ). Note that through the identification = i? i _ 1+2 «-i with i = 1, . . . , 2 9_1 , 
then q G\ = cos(0j — 0^), and using this one may stop this recursion at s = 1. 

4. Hyperspherical harmonics in polyspherical coordinates 

The general basis functions that one gets by putting coordinates on the d- dimensional 
hypersphere S d , can be specified as solutions to the angular part of Laplace's equation on 
H d+1 . These correspond to separated solutions of Laplace's equation, using the Laplace- 
Beltrami operator on the hypersphere S d . The following numbers are associated with 
each cell m G Z, I, l a , L G N (see figured] in section I3T21 . The number of vertices above 
each branching node l a and L are represented by S a and Sp respectively. The numbers 
m,l,l a ,lp label representations for the corresponding rotation subgroup in the subgroup 
chain, i.e., angular momentum type quantum numbers (see Izmest'ev et al. (2001) |16j). 

The following eigenfunctions are generated at each branching node for normalized 
hyperspherical harmonics in polyspherical coordinates using the method of trees (please 
refer to figure [1]) : 

• Type a: 

vU0 Q ) = ^=e im ^; mGZ, 0g[O,2tt). (37) 

V 27T 

• Type b: 

K h (e b ) = K> a (sm9 b ) l PP r [ a > a \cos9 b ); 

s 

n = l-lp, a = lp + Yi ^Nj, 6 b e[0,Tc]. 
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• Type b': 

< ?Q (^) = ^' /3 (cos^)^P^(sin^); I 

9 } ( 39 ) 

n = I - l a , /3 = l a + n G N , ft , G [-§ , f ] . J 

• Type c: 

*SU(*c) = 2 (Q+/3)/2+1 K' /3 (sin^)^(cos^)^Pf ' Q) (cos2^) ; ^ 

1 5 5 r ( 4 °) 

n = -(/-/ a -/ /3 ), « = / a + ^, = lp + ^-, n G N , tf c G [0, f ] . J 

We refer to the quantum number n G N in (1401) as the surrogate quantum number 
to I G N . The Jacobi polynomials P^ a '^ : C — > C can be defined in terms of the 
terminating Gauss hypergeometric series as follows (Olver et al. (2010) [2"Tl (18.5.7)]) 

p(<*,P)(~) - (" + *)» F f -n,n + « + /3+l 1- 
n\ \ a +l 2 

for n G N , and a, (3 > — 1. The normalization constants iV"'^ G (0, oo) are given by 

N a,p _ l (2n + a + (3 + l)T(n + a + (3 + l)n\ ~ 
n ~V 2 Q+/m r(n + a + l)i> + /3 + l) ' 

such that 



P^\x)Pi a ' l3 \x)(l - x) a (l + xfdx 



where a,/3 > -1 (see (7.391.1) in Gradshteyn & Ryzhik (2007) [13]). 

Notice that the eigenfunctions for cells of type b and b' can be expressed in terms 
of Gegenbauer polynomials using (see (6.4.9) in Andrews, Askey & Roy (1999) [2]) 

Cl{x) = j^^P r { r 1/2 -- l/2 \x), (41) 
where v > —1/2. Therefore we can re- write (1381) as 



^nlAVh ) = -W K—n Sint)(, "C„ ' COSb 1 !, , 42 

n ' L p y b > T(a + 1) y 2 2a+1 (2a + n)! ; v ; 

and (139]) as 



(/3 + l)y 2 2/m (2/3 + n) 
Note that even though a, (3 are not necessarily integers, 2a, 2(3 G N c 
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4-1. Examples oj hyper spherical harmonics in polyspherical coordinate systems 

With the simplest example, polar coordinates (see figure EK, Transformation ( fl9l) ). the 
normalized harmonics are given by 

p im<f> 

Y m {<j>) - 



/2tt 

In d = 3 there are two possible polyspherical coordinate systems, that of type ba 
(12Tp and b'a flU. In both cases there are two branching nodes. First we treat the type 
ba coordinate system (see figure [2b). Using (138]) we see that a = m, n = I —m, lp = m, 
and Sp = since there are no vertices above the branching node m. Through reduction 
and multiplication by the eigenf unction, the normalized spherical harmonics are 



Y Um {9,4>) = {-IT J ** 1 [) + ^j PTCoos 9)e**, (44) 

where : (—1,1) — > C is the Ferrers function of the first kind (associated Legendre 
function of the first kind on the cut) through [2T], (14.3.1)], defined as 

Notice we have used (cf. (8.936.2) in Gradshteyn & Ryzhik (2007)) [13]) 



( 1 \ m {~\ 2\—m/2 

C "™ /2(I) = (2m I 1)H Pr(x) ' (46) 

with (I42p to reduce the Gegenbauer polynomial to a Ferrers function of the first kind. 
The functions Y lm are called standard spherical harmonics. 

Now we treat the type b'a coordinate system (see figure [2b). Using f[39|) we see that 
a = m, n = I — m, l a = m, and S a = since there are no vertices above the branching 
node m. Through reduction and multiplication by the eigenf unction, the normalized 
spherical harmonics are 



YUOA) = (-lr ^^^^j PHsin^, (47) 

where we have used (]43p and flUJ) . The range of 6* in (jUJ) and fT47|) is taken to be 
9 E [0,7r] or G [— 7r/2,7r/2] (for nodes of type b and b') respectively. 

In d = 4 there are five possible polyspherical coordinate systems, those of type 
b 2 a ([25D, bb'a ([27]), b'ba ([29]), b /2 a (EI]), and ca 2 (ESJ). The normalized hyperspherical 
harmonics in type b 2 a coordinates (see figure [3^,) are 




(2l 2 + l)(l 1 + l)(l 1 -l 2 )\(l 2 -m)l 



2(/ 1 + / 2 + l)!(/ 2 + m)! 
x (sm9 1 ) l *C^l(cos9 1 )Pr 2 (cos9 2 )e im t (48 
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in type bb'a coordinates (see figure [3b) 



Y (9 9 6 ) = t i y n ^ U /(2Z a + l)(Ii + l)(Ii - f a )!(f 2 - m)! 



7r V 2(/ 1 + / 2 + l)!(/ 2 + m) 



x (sin fl^C^ (cos ^)P^(sin 9 2 )e m \ 



in type b'ba coordinates (see figure Eb) 



v (a a LjTW!! / (2/ 2 + l)(/ 1 + l)(/ 1 -/ 2 )!(/ 2 -m)! 

^A^fi, f a , 0) - - y 2(/ 1 + / 2 + l)!(/ 2 + m)! 

x (cos^) i2 CK(sm^)P^(cos^ 2 )e^, 
in type b' 2 a coordinates (see figure Eli) 



v (a a / (2/ 2 + l)(/ 1 + l)(/i-/ 2 )!(/ 2 -m)! 

^' mlSl,S,,0J ~ * V 2(/ 1 + / 2 + l)!(/ 2 + m)! 

x (cos 00^^ (sin ^)P^(sm 2 )e^, 

and in type ca 2 (see figure [3k) 
g 

^,m 1 ,m 2 (^,01,0 2 ) = - 



7T 





■ |™ll + |™ 2 |_ 




Imxl 


- Kl). 




2 W- 


\m x \ + |m 2 |) 




- I™il 


-Kl)]! 



x (sin^)l^l(cos^) |mil ^ l !;^p^ 2|)/2 (cos2^), (49) 

with the restriction to the parameter space given by \{l — \mi\ — |^ 2 |) G N . Using 
the surrogate quantum number n G N , defined in terms of / G N such that 
2n = I — Imxl — |m 2 |, we can more conveniently express the normalized hyperspherical 
harmonics in type ca 2 coordinates, namely 



V (trh ,_ e^+^) / (2w+|m 1 | + |m 2 | + l)(n+|m 1 | + |m 2 |)!n! 
^, mi , m2 W0i,0 2 J ^ y 2(n+|m 1 |)!(n+|m 2 |)! 

x (sin 0) |ma 1 (cos 0) |mi 1 P^™ 2 1 ' |mi l} (cos 20) . 

For arbitrary dimensions, we can use standard polyspherical coordinates (1351) 
to construct normalized hyperspherical harmonics. The polyspherical harmonics 
corresponding to this coordinate system are basis functions for the irreducible 
representations of 0(d) (see Vilenkin (1968) [25]). In terms of these coordinates, the 
normalized hyperspherical harmonics are 

im<j> d—2 

Y l K {x) = -^Y[ej(l j ,l j+1 ;6 j ) t (50) 

3=1 

where x G S '" 1 , K = {l 2 , l 3 , . . . , l d ^}, 

l = k>h>h>---> ld-3 > h-2 = £ > h-x = \m\ > 0, (51) 
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d-j-i 



{2l j + d-j-l){l j -l j+1 )\ 



2l j+1 + d-j -1 y wilj + lj+x + d-j -2)\ 



x (sm^^C^^-^Ccos^)- (52) 

The computation of (j52j) is a straightforward consequence of (]38p . and doing the proper 
node counting for Sg, in the tree depicted in figure HI The normalized spherical 
harmonics (PHI) are the simplest example of these harmonics. 

In generalized Hopf coordinates of type V 2 p coordinates with p > 1 (see figure 
EJ), the normalized hyperspherical harmonics : S 2 ^ 1 — > C can be given more 
conveniently using the surrogate quantum numbers n q G N which are connected to 
L G N through the relation L = l a + l B + 2n„ such that 1 < g < 2 P_1 — 1. The 
orthonormal polyspherical harmonics, which can be generated using the method of trees 
(see section I3TT1) . are 

r,f (8) = ? — - nnn , V *•) • (53) 

V27T i=1 s= i V l <xih J 

where K = {/ 2 , . . . , Lp-i,, m 1; . . . , m 2 p-i} and : Ng x [0, ~] — >■ R is defined by 



n q ^\ /(2n 9 + a + /3 + l)(n 9 + « + /?)! n g ! 



i a ,Z, ' V ' V (n g + a)!(n 9 + 0)! 



a) '(5 / y Vq 1 V<? 

x (cos # q ) la (sin ^P^' a) (cos(2tf )) 



9/ ^""gy - 1 n g \^"\.-""qJJi 

with g = 2 3 '" 1 + s - 1, l a = l 2q , l B = l 2q+1 , a = l 2q - 1 + 2 p - 2 - Llo M, an d 
/3 = ^g+i — 1 + 2 p-2 ~ L lo s 2 9J _ Note that we use the identification 

(1 2 p-i, . . . ,/ 2 p_0 = flmj , • • • , |m 2 p-i|). 

5. Jacobi polynomial and limiting expansions for the Euler kernel 

In this section we derive Jacobi, Gegenbauer and Chebyshev polynomial of the first 
kind series expansions of the Euler kernel (z — x)~ u . These series expansions are used 
to obtain azimuthal Fourier and hyperspherical harmonic expansions of a fundamental 
solution of the polyharmonic equation in Euclidean space. 

Theorem 5.1. Let a, (3 > — 1, x,z,v G C, with z G C \ (— oo, 1] on any ellipse with 
foci at ±1 and x in the interior of that ellipse. Then 

1 _ (z - l) a+1 - u (z + lf +1 ' v 
{z-x) v ~ 2 a+f3+1 ~ u 

(<* + g + 2n + !) r (« + g + " + f^p(^) f x ) f 54) 

£^ r(a + l+n)r(/3+l+n) ^n+^-l l J n I )■ K ) 
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If one substitutes z — (1 + p 2 )/(2p) in (|5"1|) with < |p| < 1, then one obtains an 
alternate expression 

1 _ (1 _ p) 2(a + l-,) (1+p) 2(^+l-,) 

(1 + p 2 - 2px) V ~ 2 2( a+P )+Z-2u pa+P+ 2~u 

^ (a + + 2n + l)r(a + + n + l)(i/) n (a+1 _„ i/?+1 _„) ( 1 + p 2 \ (0(j8) 



X 

n=0 



Proof. Consider the generating function for Gegenbauer polynomials (see for instance, 
(18.12.4) in Olver et al. (2010)) 

(1 + p -2px) ^ ^ (i/+ 5 ) n 

where we have expressed the Gegenbauer polynomial as a symmetric Jacobi polynomial 
using ( 14T|) . Now we utilize the connection relation for Jacobi polynomials (see for 
instance Ismail (2005) [HI p. 256]), namely 



>(7,<5) 
n 

k=0 



where 7, 5 > —1 and 

(7 + k + l) n - k (n + 7 + S + l) k T(a + (3 + k + l) 



Cn,fc(7,^;«,/3) 



[n-k)\T(a + f3 + 2k + l 



— n + k,n + k + j + 5 + l,a + k + 1 
j + k + l,a + /3 + 2k + 2 



The generalized hypergeometric function 3 F 2 : C 3 x (C\— N ) 2 x {z e C : \z\ < 1} — > C 
can be defined as follows 

ai,a 2 ,a 3 \ _ \^ (Ql)n( a 2)n( Q 3)n ^ 

3 " 2 ' 61,6, ' ,Z ) : ~h (PJMn nV (56) 
Utilizing (|56|) in (155]) . reversing the order of the summations, and shifting the n index 
yields after some manipulations 

1 P fc r(2z/ + 2£Qi> + (3 + fc + 1) (a><3) 

(l + p 2 -2px) ,/ 2 2l/ " 1 r(i/) ^T(z/ + A; + |)r(« + /3 + 2A;+l) fc W 

p n (2z/ + 2A;) n / -n,n + 2v + 2A;, a + fc + 1 



x 



n\ ~ \ v + k + ±a + (3 + 2k + 2 

n=0 \ z 

By substituting the definition of the 3 F 2 generalized hypergeometric function ( 1561) in the 
sum over n and reversing the orders of the resulting two summations, it follows that 

<A p n (2v + 2k) n ( -n, n + 2v + 2k, a + k + 1 



n=0 \ 2' 



v + k + ±a + (3 + 2k + 2 
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1 ( u + k,a + k + l -Ap 

72FA _ , _ , OI _, ;- -n • (57) 



"(l_p) 2 -+ 2fci 1 ^ a + / 3 + 2A; + 2 ' (1 - p) 

If we apply the right-hand side of (jSTj) to the definition of the Jacobi function of the 
second kind Q { ^ p) : C \ (-00, 1] C (cf. [TQl (10.8.18)]), namely 

2 Q+ ^r(a + 7 + l)r(/3 + 7 + l) / 7 + l,« + 7 + l 2 \ 

r(a + /3 + 2 7 + 2)(^-l) a+7+1 (^ + l) /32 1 \ « + /3 + 2 7 + 2 ' l-z) ' 1 j 

where a + 7, /3 + 7 ^ — N, noting Theorem 12.7.3 (Expansion of an analytic function in 
terms of orthogonal polynomials) in Szego (1959) we obtain the desired result. ■ 

Lemma 5.2. Let n G N , p G C \ {--, -|, — |, ...}, v G C \ -N , z 6 C \ (-00, 1]. 
o( ,-, + i/2,,-, + i/2) r , _ 2^ +1/2 r (p + n + 1/2) e ^ + V2) 1/2 

vn+^-i w r(i/ + n)(/ - i) (M_I/)/2+1/4 

Proof. Comparison of the definition of Jacobi function of the second kind fl58l) . with the 
definition of the associated Legendre function of the second kind Q„ : C \ (—00, 1] — > C 
entry 24, p. 161]) 



, e^T{v + l)T{v + p+l){z + iyl 2 ( u + l^ + p + 1 
^ [Z) - T(2u + 2)(z-lf 2+a+1 2l { 2 + 2u 'l-z 

and elsewhere in z by analytic continuation, completes the proof. ■ 

Corollary 5.3. Let v G C \ -N , with p G (-1/2, 00) \ {0}, and z G C \ (-00, 1] on 
any ellipse with foci at ±1 u>rf/i x in the interior of that ellipse. Then 

2^+l/2p/ \ i7r(Al-i/+l/2) 



If one substitutes z — (1 + p 2 )/(2p) in ( 159]) with < |p| < 1, then one obtains an 
alternate expression with x G [—1,1], namely 

1 _ r{p)e in{ ^ u+1/2) 

(1 + p 2 - 2px) v ~ 0Fr(z/)p^ +1/2 (l - p 2 )"-"- 1 ' 2 

00 / 1 i 2 ■ 



x 

n=0 



Proof. Let a = (3 = p — 1/2 in Theorem 15.11 use Lemma [5.21 and the definition of the 
Gegenbauer polynomial in terms of a symmetric Jacobi polynomial (I4ip . ■ 
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Corollary 5.4. Let d G {3, 4, 5, . . .}, u G C \ {0, 2, 4, . . .}, x, x G R d w«£/i r = ||x|| , 



r 



x — x 



x'||, and cos 7 = (x, x)/(rr') . T/ien 

e w(H-d-l)/2p ^2 _ r 2 x(^+d-l)/2 



00 / 2 , /2' 



00 

x J] (2n + d - 2) ^74/2 ) CfW), (60) 

n=0 \ 

where r< = mm {r, r'} . 

> max L ' J 



Proof. Map v 1— )■ —v/2, and then directly substitute 2; = (r 2 + r /2 )/(2rr'), x = COS7 in 
(1591) . using ||x — x'|| = V2rr' y '2 — x, x, x' G R d for <i > 3. ■ 



Corollary 5.5. Lei 1/ 6 C \ — N 0; p G (—1/2, 00) \ {0}, and x, z G C swc/i t/iat 
z G C \ (—00, 1] lie on any ellipse with foci at ±1 with x in the interior of that ellipse. 
Then 

1 \[2 e"* 71 '- 1 '" 1 / 2 - 1 00 

V^W = v^m(z 2 - ir /2 ~ 1/4 g e ^ (x)g « : v 2(z) - (61) 

If one substitutes z = (1 + p 2 )/(2p) in f lBTl) with < |p| < 1, then one obtains an 
alternate expression with x G [—1, 1], namely 

{l + p 2 -2pxY ^/WpT^il-p 2 )^ 



u-l/2 S ^2 e n T n{x)Q n -i/ 2 ( 2o ) 
n=0 ^ r / 



Proof. Take the limit as p — > on the right-hand side of ( 15T?|) and use 

Um^C^^eJ,,^) (62) 

(see for instance (6.4.13) in Andrews, Askey & Roy (1999) |2j), where T n : C — > C 
is the Chebyshev polynomial of the first kind defined as (see section 5.7.2 in Magnus, 
Oberhettinger & Soni (1966) 

T n (x) := 2 F 1 



-n,n 1 - x 



Note that (EE]) is given in Cohl & Dominici (2011) |5l (3.10)], so (J53H and (j5S) represent 
generalizations of that formula. 
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6. Addition theorems in Vilenkin's polyspherical coordinates 

As mentioned in the introduction, by comparing azimuthal Fourier expansions for 
fundamental solutions of the polyharmonic equation on R d for d > 3, with eigenfunction 
expansions for these fundamental solutions in rotationally-invariant coordinate systems, 
one can construct multi-summation and multi-integration addition theorems. 

We have already constructed some of these addition theorems on R 3 (see Cohl 
et al. (2000) [8]; Cohl et al. (2001) (7J). In Cohl et al. (2000) [8] addition theorems 
are derived in this fashion for fundamental solutions of the Laplacian on R 3 , namely 
in cylindrical, oblate spheroidal, prolate spheroidal, parabolic, bispherical and toroidal 
coordinate systems. In Cohl et al. (2001) [7J an addition theorem is derived from a 
fundamental solution of Laplace's equation in type ba (standard) spherical coordinates 
(we generalize that addition theorem below 11731)). In circular cylindrical and rotational 
parabolic coordinates, the derived addition theorems are given in terms of definite 
integrals. In oblate spheroidal, prolate spheroidal, bispherical, toroidal, and spherical 
coordinates, the derived addition theorems are given in terms of sums. As is indicated 
by the above presentation, one may construct addition theorems in this manner in 
any rotationally-invariant coordinate system which yields solutions through separation 
of variables for the polyharmonic equation. In a similar setting, addition theorems 
may be generated for other inhomogeneous linear partial differential equations, such 
as for the Helmholtz, wave and heat equations in arbitrary dimensions. Furthermore 
an extension of this concept is possible when working with linear partial differential 
operators on Riemannian manifolds, such as for the Laplace-Beltrami operator (see Cohl 
& Kalnins (2012) [61 section 5.1]). Once a Fourier expansion for a fundamental solution is 
obtained for a partial differential operator on a Riemannian manifold, one must construct 
eigenfunction expansions for a fundamental solution corresponding to that operator and 
identify those nested multi-summation and multi-integration eigenfunction expansions 
which correspond to the Fourier coefficient for that operator. 

Suppose one adopts for R d a transformation using standard polyspherical 
coordinates (135]) . New addition theorems can then be generated by re-arranging the 
relevant sums over the hyperspherical harmonics, for instance 

^1 max ^Imax ^lmax ^lmax ^lmax ^lmax 

EE = E E E -"EE- («) 

Z 1= K ™=-hmtn. l d-2=\ m \ l d-3=h-2 h=h h= l 2 

In ( )63|) . one must sum over all allowed quantum numbers {1{\UK. This is accomplished 
by taking the limit on both sides of ( )63|) as l\ max — > oo. The forward sum on the left-hand 
side of (163j) is now to be interpreted in a reverse fashion as in the right-hand side. Take 
for instance the limit ^i max — > oo in the final sum, i.e., {/ x G N : Z x > l 2 }. Once l x has 
been summed over and eliminated from the sum, then one must sum over the remaining 
quantum numbers. In this limit, one then sums over l 2 in the set {l 2 G N : l 2 > Is}, 
eliminating l 2 from the sum. One continues this process until we have eliminated the 
contribution from all quantum numbers except m and finally we are left with a series 
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which only depends on m. Since we have already identified the azimuthal Fourier 
coefficients in section [5] (the terms corresponding to m with the other quantum numbers 
completely summed over), then by equating these two expressions, we can construct 
multi-index addition theorems in any polyspherical coordinate system. In the discussion 
below, we give explicit examples of this procedure. 

6.1. Power-law addition theorem on R d for d > 3 in standard polyspherical coordinates 
In standard polyspherical coordinates (135 p we have the following addition theorem. 

Theorem 6.1. Let v E C\{2m, 2m+2, 2m+4, . . .}, m E Z, B i E [0, ir], for 1 < i < d-2, 
r, r E [0, oo) with d > 3. Then 

^ ( -T/ 2 )/2 (x) = V2e^ d ~ 2 ^ d ~^ 2 ^rr'nsin^sin^' 

x (V - 1) ~ {u+m l r > z r <) 

X [X 1} (rr'f- 1)/2 

x / ; r; r, cost^ 9 r ; cos 6*, 9 

^ (/ d _ 2 + m)! ld ~ 2K d ~ 2J ld ~ 2K d ~ 2> 



oc 

X 

'd-3='d-2 



©d-3 (^Z-3; ^d-2i ^d-3) @d-3 (^d-3; ^d-2j ^-3) 



00 

X 



E 2 fe/3^2) (Za./aJ^) 



00 / 2 . 

,(i-^-d)/2 ( r + r 

!+(d-3)/2 I 2rr ' 



x £ 9? (/ l5 / 2 ; AO 9? Z 2 ; 6»i) , (64) 



where 



d-2 t-1 

2 , /2 



r 



+ r' — 2rr'^^ cos ^cos ^ J^J sin ^ sin 9j 

i=l j= 
d-2 

2rr' J^J sin ^sin Q\ 



X = ^ • (65) 



i=i 



Proof. If we adopt standard polyspherical coordinates fl35|) (see figure H]), then one 
can obtain an eigenfunction expansion for a power-law fundamental solution of the 
polyharmonic equation in standard polyspherical coordinates using the Gegenbauer 
expansion (|60|) with the addition theorem for hyperspherical harmonics ([IT), (150)) . The 
Fourier expansion for a power-law fundamental solution of the polyharmonic equation 
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is obtained from the algebraic expression for a power-law fundamental solution of the 
polyharmonic equation (cf. ([7])) and (1611) . namely 

!— in(u+l)/2 ( d-2 \ v / 2 



x — X 



x £ e^)Q^+r(x), (66) 

m=— oo 

where 



r 



d-2 i-1 

2 + r' 2 — 2rr'^^ cos #jCos 9[ J^J sin 9j sin ^- 



X = =^ ^ ■ («) 



2rr' J^J sin ^sin Q\ 



8=1 



The Gegenbauer expansion yields 

4e i W (i/+d-i)/2 7r (d-i)/2 ^2 _ r 2^(^+<i-i)/2 



x — X 



A (d-\)/2 



r( 2) (rr')' 



x 

Now if we expand the product of polyspherical harmonics in (1681) with (1501) . ( |63l) . we 
obtain 



x - x 'ir = Yl c ° s ( m (^ - — v(_u\ 

m=0 \ 2) 



-3)/2 e ^-D/2 em ( r 2 _ r 2 ) (^-l)/2 



(d-l)/2 



00 



x ^ + ^ - m)! pg_, ( c«^) pil (cose,) 

Id--?— m 



X 

ld-Z=ld-2 



@d-3 ^2-2) ©d-3 (^d-3> ^2-2> ^d-s) 



oo 

^2='3 

oo 2 , /2 

< £ e? (h, i 2 - e x ) ei {h, i 2 - 9[) Q?-\;% {-^r-) , (69) 

/ 1= / 2 



X;eg(z 2 ,z 3 ;0 2 )eg (/ 2 ,/ 3 ^2) 



where 0^(L, L +1 ;0), for 1 < j < d — 2 is defined in fl52|) . By comparing the Fourier 
coefficients of fl69|) with (166]) . we complete the proof of this theorem. ■ 
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This is just one example of a derived multi-summation addition theorem for arbitrary 
dimensions. There are an unlimited number of such straightforward examples to 
generate. In the next section we derive another example which is valid on R d where d 
is given by a power of two, i.e., generalized Hopf coordinates (|36|) . 

2 q 

6.2. Power-law addition theorem on R for q > 2 in generalized Hopf coordinates 

Theorem 6.2. Let v G C \ {2m, 2m + 2, 2m + 4, . . .}, m 1 G Z, r,r' G [0,oo), 4 G [0, f ] 
tcita 1 < i < 2 9 " 1 - 1, 0i G [0, 2tt) suc/i that l<i< 2 9_1 . Then 

'q-i \ -»/ 2 

J cos fig-i cos &2 



V mi _l/ 2 U; 2 (v+l)/2^2 _ 1 ^+l)/4 



2 N (^+2 9 -l)/2 



/ 



oo 

^ e ma cos(m 2 (0 2 - 0' 2 )) • • • e m cos(m 29 -i(0 29 -i - 0' 2 ,_i)) 



2 4 

m 2 =0 m 29 _i=0 



oo 

X 



2 |m 29 -i_J, |m 2 ,-i| y |m 2 ,-i_J, |m 2 ,-i| 



V-i_i =0 



— „ V m, , m, / V m, , m 9 

n 2<3 _2=0 v 1 11 1 ' v 1 11 1 



oo 



V- 2 -l 



=0 



oo / x/ -, / 2 , 

Ui ■Ao (i - v - 2 " )/2 f r +r 

J J ' 1 y V 2n 1 +Z 2 +Z 3 + (2 9 -3)/2 I 2 r/ 



n 1= 

9-1 

r 2 + r' 2 — 2rr'cos7 + 2rr'cos(0 1 — </)[) JJcostfjj-i cos$^-i 

X = i=i • (71) 

2rr' cos ?9 2 j-i cos 
i=i 

Proof. If we adopt generalized Hopf coordinates fl36l) (see figure [5]), then we can 
use the corresponding harmonics (1531) in combination with the addition theorem for 
hyperspherical harmonics (Q. We compare the Gegenbauer expansion for powers of the 
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distance (IBP]) with the Fourier expansion 

/ q -i \ "I 2 

^2ie inv/2 2rr' JJ cos tf^-i cos 

>uv V j=l 



X — X 



v^r (=*) ( X 2 - i)-^ 4 

oo 

E j^^Q^-i/fix), (72) 



x 

m-i =— oo 



where x > 1 is given by (ITT]) . Notice that x is independent of 4>i — 4>i- By using the 
Gegenbauer expansion (160]) and inserting the appropriate Gegenbauer polynomial using 
the addition theorem for hyperspherical harmonics ([T|), we obtain 



2 q 1 iwu/2/ 2 2\0+2 ? -l)/2 



x — X = - 



?r(^)K) <! "- 1,/2 

x E CS-w, (^) (S) W)- (73) 

By expanding the product of polyspherical harmonics in (173"]) with fl53|) expressed in 
terms of the surrogate quantum number n and reversing the order of the sums, we obtain 
a multi-summation expression for the power of the Euclidean distance between two 
points in generalized Hopf coordinates. Through comparison of the resulting equation 
with the m 1 Fourier coefficients of (172]) . we derive (170]) . a mult i- summation addition 
theorem for the associated Legendre function of the second kind with argument x 
(cf. (J7JJ)). ■ 

6.3. Power-law addition theorems on R 3 

In d = 3 there are two ways to construct polyspherical coordinates, with trees of type 
b'a and ba. 

6.3.1. Type ba coordinates 

Corollary 6.3. Let v e C \ {2m, 2m + 2, 2m + 4, . . .}, m e N 0; 9,9' e [0,tt], 
r, r G [0, oo). T/jen 

/ 2 2 \ 0+2)/2 

^V2 )/2 (X) = ^2-^ 3 )/ 2 (sin^sin^)^ /2 (x 2 - 1)^ +1)/4 ( ! ^) 

(f-m)l ( , +2)/2 /r^+r' 2 
(Z + m)! ' I 2rr 



x E< 2 ' + Dfe^T^^ 2 " 2 ?T^)??{^% (74) 



w/iere 

2 , /2 



r + r' — 2rr' cos 9 cos 

X = /•/)■/)/ • ( 75 ) 

Irr sin t/ sin 9 
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Proof. Taking d = 3 in (1641) completes the proof. ■ 

Equation (174)) is a generalization of one of the main results of Cohl et al. (2001) [7J. This 
can be observed if you substitute v — —1 in ( 1741) (this corresponds to a fundamental 
solution of Laplace's equation on R 3 ), then the associated Legendre function of the 
second kind on the right-hand side reduces to an elementary function through (TJ 
(8.6.11)], producing 




6.3.2. Type ba coordinates 

Corollary 6.4. Let v G C \ {2m, 2m + 2, 2m + 4, . . .}, m G Z, 9,9' G [-f,f], 
r, r G [0, oo). T/ien 

/ 2 2 \ (H-2)/2 

Q-J-$ ,2 M = ! v^2-<"«»/ 2 (coBeco S er" /2 (x 2 - ir ( " +1)/1 ( ri ^ ri ) 



X £< 2 < + ^ra^" 2 Pr(3inO)Pr(sinO, (77) 



i=|m| 

where 

r 2 + r' 2 — 2rr sin 6* sin 6*' 
^ 2rr' cos 6* cos 9' 

Proof. Using the method of trees, the map from type ba to type b'a coordinates 
is given by 9, 9' i— > 9 — |. Under these transformations cos 9, cos 9' \— > sin 9, sin 0' and 
sin 0, sin 0' i— > — cos 0, — cos 9' , as applied to ( 1741) . This completes the proof. ■ 

fr^. Power-law addition theorems on R 4 

In d = 4 there are five ways to construct polyspherical coordinates, with trees of type 
b 2 a (see figure Et), bb'a (see figure [3b) , b'ba (see figure Eb), b /2 a (see figure Eli), ca 2 
(see figure Efe)- 

6.4..I. Type b 2 a coordinates 

Corollary 6.5. Let v G C \ {2m, 2m + 2, 2m + 4, . . .}, m G Z, r, r' G [0, 00), 

0i,0i,0 2 ,0a e [0,4 Then 

* / 2 2 n 0+3)/2 

00 02*0/07 , lW7 l\2/ 



x > — (sm0 1 sm0 1 2 P, cos 2 P 2 2 (cos0 2 ) 

(l 2 + my. 

h=\m\ 
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H 11 + l)( h -l 2 y _ {v+3)/2 ( r 2 + r > r 
(h + l 2 + l)\ ^ +1 / 2 I 2rr' 



x (] i / i m ^1+1/2 w ^^(cosgiJ^-^Ccosg!), (7? 



where 

r 2 + r' 2 — 2rr' cos 9 X cos 9[ — 2rr sin 6 1 ! sin 9[ cos # 2 cos $ 2 
2rr' sin sin 9[ sin 6* 2 sin 9' 2 

Proof. Taking d = 4 in f !64p completes the proof. ■ 

If you substitute v = —2 (a fundamental solution for the Laplacian on R 4 ) in (I7HI) then 
the Legendre functions of the second kind reduce to elementary functions through [TJ 
(8.6.10-11)], and one obtains the following 

, — = 2 sin #1 sin #1 sin 9 2 sin #9 

x E ^ 2 " (5ingl Bin^)' 2 P^(cosg 2 )P^(cosg 2 ) 

i 2 =|m| 

6.4-2. Type bb'a coordinates In type bb'a coordinates, this is the addition theorem 
derived using the methods described above (see figure [3b) 

Corollary 6.6. Let v G C\{2m, 2m+2, 2m+4, . . .}, m G Z, r, r G [0, 00), 9 x ,9'i G [0, tt], 
2 ,0 2 G[-f,§]. Then 

, / 2 2 \ 0+3)/2 

= ^ (^) (X 2 - l)- ( " +I,/4 (sin« 1 sin< ) ' 1 cos0 2 cos^)-'- /2 

~ 2 2 ' 2 (2/ 2 + l)(/ 2 !) 2 (/ 2 -m)! . w 2T3m , . fl , pm/ . 

x V 77— ^ sm^sm^ 2 P, sm# 2 P« 2 sm# 2 

( 2 =|m 

(ji + l)(/j - ^)! (y+3 )/a / r 2 + r' 2 ^\ 



X 

w/iere 

r 2 + r' 2 — 2rr' cos 9 1 cos 6*i — 2rr' sin 9 1 sin #1 sin 9 2 sin # 2 

X = n I ■ n ■ nl n 7? ■ ( 81 ) 



2rr sin 9i sin 9 1 cos 6* 2 cos 9 



2 



Proof. Using the method of trees, the map from type b 2 a to type bb'a coordinates 
leaves 9 1 invariant and 9 2 ,9' 2 i— > 9 2 — f,0 2 — §. Under these transformations 
cos 9 2 , cos 9' 2 i — ^ sin ^2, sin ^ 2 and sin #2, sin # 2 ^ — cos# 2 , — cos ^2, as applied to ( 180]) . 
This completes the proof. ■ 
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6.4-3. Type b'ba coordinates In type b'ba coordinates, this is the addition theorem 
derived using the methods described above (see figure [3b) 

Corollary 6.7. Let v G C \ {2m, 2m + 2, 2m + 4, . . .}, m G Z, r, r G [0, oo), 
0i,0i G [-f,|], 2 ,0 2 e [0,4 Then 

, / 2 2 \ 0+3)/2 

Q^TW = ^ (^) (X 2 - ^-^(cos^cos^sm^sin^)-^ 

o2«o /Ol |\2/ 



^2 2 (2/ 2 + l)(/ 2 !) {h-m)\, Q f\h-om r n \pm/ fl /- 

X > — COS 0! COS 0! 2 P, (COS 2 )P Z COS 2 

(< 2 + mj! 



z 2 =H 

V"^ (*1 + ~ kV- n -(u+3)/2 I T 2 +r \ , +1 s^a+l / - fl /x /o 9 n 

X 2_, (/ 1 + / 2 + i)i ^ii+i/2 ^ 2rr' J C h-h( sme V C ii-i2( sm6 U> ( 82 > 

where 

r 2 + r' 2 — 2rr' sin 6 1 sin 0^ — 2rr' cos 0! cos 0^ cos 2 cos 2 . . 
2rr cos X cos 9 1 sin 2 sin 2 

Proof. Using the method of trees, the map from type b 2 a to type b'ba coordinates 
leaves 2 invariant and 0i,0i h- >■ 6 l — |,0 X — |. Under these transformations 
cos X , cos 0'x i — y sm9 1 ,sm9[ and sin0 1 ,sin0' 1 h-> — cos0 l7 — cos0i, as applied to ( IHUj) . 
This completes the proof. ■ 

6.4-4- Type b /2 a coordinates In type b /2 a coordinates, this is the addition theorem 
derived using the methods described above (see figure [3bl) 

Corollary 6.8. Let v G C \ {2m, 2m + 2, 2m + 4, . . m <E Z, r,r' e [0, oo), 
0i,0'i,02,0 2 e [-§,§]. Then 

Qm-lft^lx) = 2 (,+l)/2 ( r> rr / r< j " 1) _(!/+1)/4 (COS0! COS 9[ COS0 2 COS0 2 )^ /2 

r>2«, /r>7 , |\2, 



(_,: 



--\m\ 



v 2^(2/ 2 + l)(/ 2 !)^(/ 2 -m)! / , , . . 

x y. 77—, r, (cos0 1 cos0 1 ) 2 P, (sm0 2 )P, 2 (sm0 2 

(t 2 + mj! 



(Zi + l)(Zi-Z 2 )!_ _ (v+3 )/2 r+r' ) 



1—^2 



X 

u>/iere 

r 2 + r' 2 — 2rr' sin X sin 9[ — 2rr' cos X cos 0' x sin 2 sin 9' 2 
2rr cos 9 l cos X cos 2 cos 2 

2 /2 

Proof. Using the method of trees, the map from type b a to type b a 
coordinates 9 1 ,9' 1 ^ 9 l - § , 9[ - § and 2 , 2 ^ 2 - f, 9' 2 - §. Under 
these transformations cos 9 1 , cos 9[ , cos 2 , cos 9' 2 H- sin X , sin 9[ , sin 2 , sin 9' 2 and 
sin 1 ,sin0 1 h> — cos0 1; — cos0' l7 — cos0 2 , — cos0 2 , as applied to flHDj) . This completes 
the proof. ■ 
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6.4-5. Type ca 2 coordinates In type ca 2 coordinates, this is the addition theorem 
derived using the methods described above 

Corollary 6.9. Let v G C \ {2m, 2m + 2, 2m + 4, . . .}, m 1 G Z, r, r G [0,oo), 
G [0,f], 2 ,0' 2 G [0,2tt). TTien 



-(i/+l)/2 / 2 2 \ (H-3)/2 

0^(x) - ^ (^) (cos W) 1 ™''-^ 



oo 

l\m 2 



x e m2 cos(m 2 (0 2 - 2 ))(sin$sin$ 

m 2 =0 
oo 

71=0 



E 



(2n + 


mi 


+ m 2 + l)(|mi 


+ m 2 + n)\n\ 


(|mi 


+ n)!(m 2 + n)! 



x f ) P^' l " lll) (co S 2^)Pi^l^l)(co S 2^), (86) 



where 



r 2 + r' 2 — 2rr' sin d sin i?' cos(0 2 — 2 ) 

^ = / q ^ ' ( 87 ) 

2rr cos # cos v 

Proof. Taking q = 2 in (!7D|) completes the proof. ■ 

If you substitute v = —2 in (186]) . then the Legendre functions of the second kind reduce 
to elementary functions through [U (8.6.10-11)], and one obtains the following 

= 2 (cos d cos ' mi ' +1 e m2 cos [m 2 (0 2 — 2 )] (sin $ sin -i?')'™ 2 

(X+VX - I)™' m 2 =0 



n=0 



E 



(|m! 


| + 


m 2 


+ n)!n! 


(|mi| 


+ n 


)!(|m 2 


+ n)! 



X 



P^I-KI)(cos2^)Pi |m2l,|mil) (cos2^). (88) 



Note that in the addition theorems (1861) and ( 1881) . that if you make the map d H- <& — |, 
then this transformation preserves the addition theorems such that rri\ m 2 . This 
transformation is equivalent to swapping the position of (f>i and 2 for the tree in figure 
Efe. 
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